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We propose a quantum algorithm for solving combinatorial search problems that uses only a se- 
quence of measurements. The algorithm is similar in spirit to quantum computation by adiabatic 
evolution, in that the goal is to remain in the ground state of a time-varying Hamiltonian. Indeed, 
we show that the running times of the two algorithms are closely related. We also show how to 
achieve the quadratic speedup for Grover's unstructured search problem with only two measure- 
ments. Finally, we discuss some similarities and differences between the adiabatic and measurement 
algorithms. 



I. INTRODUCTION 

In the conventional circuit model of quantum compu- 
tation, a program for a quantum computer consists of a 
discrete sequence of unitary gates chosen from a fixed set. 
The memory of the quantum computer is a collection of 
qubits initially prepared in some definite state. After a 
sequence of unitary gates is applied, the qubits are mea- 
sured in the computational basis to give the result of the 
computation, a classical bit string. 

This description of a quantum computer has been used 
to formulate quantum algorithms that outperform clas- 
sical methods, notably Shor's factoring algorithm Q and 
Grover's algorithm for unstructured search [||. Subse- 
quent development of quantum algorithms has focused 
primarily on variations of the techniques introduced by 
Shor and Grover. One way to motivate new algorithmic 
ideas is to consider alternative (but in general, equiva- 
lent) descriptions of the way a quantum computer oper- 
ates. For example, the technique of quantum computa- 
tion by adiabatic evolution || is most easily described by 
a quantum computer that evolves continuously according 
to a time- varying Hamiltonian. 

Another model of quantum computation allows mea- 
surement at intermediate stages. Indeed, recent work has 
shown that measurement alone is universal for quantum 
computation: one can efficiently implement a universal 
set of quantum gates usingonly measurements (and clas- 
sical processing) j|, || [| [7). In this paper, we describe 
an algorithm for solving combinatorial search problems 
that consists only of a sequence of measurements. Using 
a straightforward variant of the quantum Zeno effect (see 
for example || ||, [To)), we show how to keep the quan- 
tum computer in the ground state of a smoothly varying 
Hamiltonian H(s). This process can be used to solve a 
computational problem by encoding the solution to the 
problem in the ground state of the final Hamiltonian. 

The organization of the paper is as follows. In Section 
II, we present the algorithm in detail and describe how 



measurement of H(s) can be performed on a digital quan- 
tum computer. In Section HI, we estimate the running 
time of the algorithm in terms of spectral properties of 
H(s). Then, in Section IV, we discuss how the algorithm 
performs on Grover's unstructured search problem and 
show that by a suitable modification, Grover's quadratic 
speedup can be achieved by the measurement algorithm. 
Finally, in Section [v|, we discuss the relationship between 
the measurement algorithm and quantum computation 
by adiabatic evolution. 



II. THE MEASUREMENT ALGORITHM 

Our algorithm is conceptually similar to quantum com- 
putation by adiabatic evolution ||, a general method 
for solving combinatorial search problems using a quan- 
tum computer. Both algorithms operate by remaining 
in the ground state of a smoothly varying Hamiltonian 
H(s) whose initial ground state is easy to construct and 
whose final ground state encodes the solution to the prob- 
lem. However, whereas adiabatic quantum computation 
uses Schrodinger evolution under H(s) to remain in the 
ground state, the present algorithm uses only measure- 
ment of H(s). 

In general, we are interested in searching for the mini- 
mum of a function h{z) that maps n-bit strings to positive 
real numbers. Many computational problems can be cast 
as minimization of such a function; for specific examples 
and their relationship to adiabatic quantum computa- 
tion, see |3| [OJ. Typically, we can restrict our attention 
to the case where the global minimum of h(z) is unique. 
Associated with this function, we can define a problem 
Hamiltonian Hp through its action on computational ba- 
sis states: 



H P \z)=h{z)\z) 



(1) 



Finding the global minimum of h(z) is equivalent to find- 
ing the ground state of Hp. If the global minimum is 
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unique, then this ground state is nondegenerate. 

To reach the ground state of Hp, we begin with the 
quantum computer prepared in the ground state of some 
other Hamiltonian Hp, the beginning Hamiltonian. Then 
we consider a one-parameter family of Hamiltonians H(s) 
that interpolates smoothly from Hb to Hp for s <E [0, 1]. 
In other words, H(0) — Hb and H(l) = Hp, and the 
intermediate H (s) is a smooth function of s. One possible 
choice is linear interpolation, 



H(s) = (1 - s)H B + sH P . 



(2) 



Now we divide the interval [0, 1] into M subintervals of 
width S = 1/M. So long as the interpolating Hamiltonian 
H(s) is smoothly varying and 5 is small, the ground state 
of H(s) will be close to the ground state of H(s + S). 
Thus, if the system is in the ground state of H (s) and 
we measure H(s + S), the post-measurement state is very 
likely to be the ground state of H(s + S). If we begin 
in the ground state of H(0) and successively measure 
H(S), H(25), H((M ~ 1)5), H (I), then the final state 
will be the ground state of H (1) with high probability, 
assuming 5 is sufficiently small. 

To complete our description of the quantum algorithm, 
we must explain how to measure the operator H(s). The 
technique we use is motivated by von Neumann's descrip- 
tion of the measurement process ||. In this description, 
measurement is performed by coupling the system of in- 
terest to an ancillary system, which we call the pointer. 
Suppose that the pointer is a one-dimensional free parti- 
cle and that the system-pointer interaction Hamiltonian 
is H(s) ® p, where p is the momentum of the particle. 
Furthermore, suppose that the mass of the particle is 
sufficiently large that we can neglect the kinetic term. 
Then the resulting evolution is 



-itH(s)®p 



^[|K(s))(K(s)|®e 



-UE a (s)p 



(3) 



where \E a (s)) are the eigenstates of H(s) with eigenval- 
ues E a (s), and we have set H — 1. Suppose we prepare the 
pointer in the state \x — 0), a narrow wave packet cen- 
tered at x = 0. Since the momentum operator generates 
translations in position, the above evolution performs the 
transformation 

\E a (s)) ®\x = 0)^ \E a (s)) ®\x = tE a (s)) . (4) 

If we can measure the position of the pointer with suf- 
ficiently high precision that all relevant spacings x a b — 
t\E a (s) — Eb(s)\ can be resolved, then measurement of 
the position of the pointer — a fixed, easy-to-measure 
observable, independent of H(s) — effects a measure- 
ment of H(s). 

Von Neumann's measurement protocol makes use of 
a continuous variable, the position of the pointer. To 
turn it into an algorithm that can be implemented on 
a fully digital quantum computer, we can approximate 
j) using r quantum bits to represent the 
. The full Hilbert space is thus a tensor 



the evolution 
pointer pi Il3 



product of a 2™-dimensional space for the system and a 
2 r -dimensional space for the pointer. We let the com- 
putational basis of the pointer, with basis states 
represent the basis of momentum eigenstates. The label 
z is an integer between and 2 r — 1, and the r bits of 
the binary representation of z specify the states of the r 
qubits. In this basis, the digital representation of p is 



E 2 " 

i=i 



(5) 



a sum of diagonal operators, each of which acts on only 
(i) 

a single qubit. Here al is the Pauli z operator on the 
jth qubit. As we will discuss in the next section, we have 
chosen to normalize p so that 



pz) 



(6) 



which gives ||p|| ~ 1. If H(s) is a sum of terms, each of 
which acts on at most k qubits, then H (s) ® p is a sum 
of terms, each of which acts on at most k + 1 qubits. As 
long as k is a fixed constant independent of the problem 
size n, such a Hamiltonian can be simulated efficiently on 
a quantum computer [ fl4"| . Expanded in the momentum 
eigenbasis, the initial state of the pointer is 



f 



2-1 



l* = °) = ^2» 



(7) 



The measurement is performed by evolving under H(s)(3 
p for a total time r. We discuss how to choose r in the 
next section. After this evolution, the position of the 
simulated pointer could be measured by measuring the 
qubits that represent it in the x basis, i.e., the Fourier 
transform of the computational basis. However, note that 
our algorithm only makes use of the post-measurement 
state of the system, not of the measured value of H(s). In 
other words, only the reduced density matrix of the sys- 
tem is relevant. Thus it is not actually necessary to per- 
form a Fourier transform before measuring the pointer, 
or even to measure the pointer at all. When the system- 
pointer evolution is finished, one can either re-prepare 
the pointer in its initial state \x = 0) or discard it and 
use a new pointer, and immediately begin the next mea- 
surement. 

As an aside, note that the von Neumann measurement 
procedure described above is identical to the well-known 
phase estimation algorithm for measuring the eigenvalues 
of a unitary operator |TJ|, [lTJ] , which can also be used to 
produce eigenvalues and eigenvectors of a Hamiltonian 
JT7[ | . This connection has been noted previously in [ p"3| |, 
and it has been pointed out that the measurement is a 
non-demolition measurement in Ea |. In the phase es- 
timation problem, we are given an eigenvector of a 
unitary operator U and asked to determine its eigenvalue 
e -»0_ rp^g algorithm uses two registers, one that initially 
stores \ip) and one that will store an approximation of 
the phase (f>. The first and last steps of the algorithm are 
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Fourier transforms on the phase register. The intervening 
step is to perform the transformation 



\ip) <8> \z) -> U z \i>) <g> \z) 



(8) 



where \z) is a computational basis state. If we take \z) 
to be a momentum eigenstate with eigenvalue z (i.e., if 
we choose a different normalization than in (Q)) and let 
U = e~ lHt , this is exactly the transformation induced by 
e -i(H®p)t_ ^hus we see -that the phase estimation algo- 
rithm for a unitary operator U is exactly von Neumann's 
prescription for measuring i In U. 



III. RUNNING TIME 

The running time of the measurement algorithm is the 
product of M, the number of measurements, and r, the 
time per measurement. Even if we assume perfect projec- 
tive measurements, the algorithm is guaranteed to keep 
the computer in the ground state of H(s) only in the 
limit M — > oo, so that 5 — l/M — > 0. Given a finite 
running time, the probability of finding the ground state 
of Hp with the last measurement will be less than 1. To 
understand the efficiency of the algorithm, we need to 
determine how long we must run as a function of n, the 
number of bits on which the function h is defined, so that 
the probability of success is not too small. In general, if 
the time required to achieve a success probability greater 
than some fixed constant (e.g., |) is poly(n), we say the 
algorithm is efficient, whereas if the running time grows 
exponentially, we say it is not. 

To determine the running time of the algorithm, we 
consider the effect of the measurement process on the 
reduced density matrix of the system. Here, we simply 
motivate the main result; for a detailed analysis, see Ap- 
pendix W. 

Let pM* denote the reduced density matrix of the sys- 
tem after the jth measurement; its matrix elements are 



Pab 



(E a (j5)\pW\E b (jS)) 



(9) 



The interaction with the digitized pointer effects the 
transformation 

\E a (s)) ® \z) - e -^ zt l r \E a {s)) ® \z) . (10) 

Starting with the pointer in the state (Q), evolving ac- 
cording to (pp|), and tracing over the pointer, the quan- 
tum operation induced on the system is 



Pab 



/ j u ac Pcd u bd 



(11) 



where the unitary transformation relating the energy 
eigenbases at s = jS and s = (j + 1)£ is 



and 



U$ = (E a ((j + l)5)\E b (j6)) 

- y 



J[E b (jS)-E a (j6)]zt/2 r 



(12) 



(13) 
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FIG. 1: The function \k{x)\ 2 for r = 4. 



Summing this geometric series, we find 



.U) 

"ab 



\K{[E b (j5)-E a {j5))t/2)Y 



where 



l^)| 2 = 



sin 2 x 



(14) 



(15) 



4 r sin 2 (a;/2 r ) ' 

This function is shown in Fig. |l| for the case r = 4. It has 
a sharp peak of unit height and width of order 1 at the 
origin, and identical peaks at integer multiples of 2 r ir. 

If the above procedure were a perfect projective mea- 
surement, then we would have K ab — whenever E a ^= 
E b - Assuming (temporarily) that this is the case, we find 



Poo +1) > 



U, 



Poo 



(16) 



with the initial condition Pqq = 1 and pfj} — otherwise. 
Perturbation theory gives 



U, 



= 1-5 



\(E a (s)m\E (s))f 



(Eo(s)~E a (s)r 



s=jS 



> 1 



+0(5 3 ) 

r(jsfs 2 



0(6 3 



(17) 
(18) 



g(jS) 2 

where 

T( S ) 2 = (E (s)\(^) 2 \E (s)) - (E ( s )\^\E (s)r (19) 



and 



;(«) = £i(s) - E (s) 



(20) 



is the energy gap between the ground and first excited 
states. If we let 



2 = 



r = max r(s) 

«6[0,X] 

q = min q(s) , 

«6[0,X] 



(21) 
(22) 
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then according to ([Tq), the probability of being in the 
ground state after the last measurement is at least 



(M) . 

Poo > 



-, M 



exp 



M 2 5 2 

r 2 



Mr 



+ 0(M~ 3 ) 
+ 0(M~ 2 ) . 



The probability of success is close to 1 provided 



r 2 

M > 

5 



(23) 
(24) 

(25) 



When if a and ifp are both sums of poly(n) terms, each 
of which acts nontrivially on at most a constant number 
of qubits, it is easy to choose an interpolation such as 
(Q) so that r is only poly(n). Thus we are mainly inter- 
ested in the behavior of g, the minimum gap between the 
ground and first excited states. We see that for the algo- 
rithm to be successful, the total number of measurements 
must be much larger than 1/g 2 . 

In fact, the simulated von Neumann procedure is not 
a perfect projective measurement. We must determine 
how long the system and pointer should interact so that 
the measurement is sufficiently good. The analysis in Ap- 
pendix [A] shows that |«oi | 2 should be bounded below 1 
by a constant for all j. In other words, to sufficiently re- 
solve the difference between the ground and first excited 
states, we must decrease the coherence between them by 
a fixed fraction per measurement. The width of the cen- 
tral peak in Fig. [j] is of order 1, so it is straightforward 
to show that to have |«(a;)| 2 less than, say, 1/2, we must 
have x > 0(1). This places a lower bound on the system- 
pointer interaction time of 



r > 



0(1) 



9 



(26) 



independent of r, the number of pointer qubits. (Note 
that the same bound also holds in the case of a continuous 
pointer with a fixed resolution length.) 

Putting these results together, we find that the mea- 
surement algorithm is successful if the total running time, 
T = Mr, satisfies 



r 2 

T>-r (measurement) . 
9 A 



(27) 



This result can be compared to the corresponding expres- 
sion for quantum computation by adiabatic evolution, 



9' 



(adiabatic) . 



(28) 



Note that the same quantity appears in the numerator 
of both expressions; in both cases, T accounts for the 
possibility of transitions to all possible excited states. 

The adiabatic and measurement algorithms have qual- 
itatively similar behavior: if the gap is exponentially 
small, neither algorithm is efficient, whereas if the gap 



is only polynomially small, both algorithms are efficient. 
However, the measurement algorithm is slightly slower: 
whereas adiabatic evolution runs in a time that grows 
like 1/g 2 , the measurement algorithm runs in a time that 
grows like 1/g 3 - To see that this comparison is fair, re- 
call that we have defined the momentum in @ so that 
~ 1, which gives ||-ff(s)|| ~ \\H(s) Alterna- 
tively, we can compare the number r\ of few-qubit uni- 
tary gates needed to simulate the two algorithms on a 
conventional quantum computer. Using the Lie product 
formula 



a A+B 



A j ra / myn 



(29) 



which is valid provided m ^> \\A\\ + ||B|| , we find n = 
0(1/ g 4 ) for adiabatic evolution and n = 0(1/ g 6 ) for the 
measurement algorithm, in agreement with the previous 
comparison. 

The argument we have used to motivate (|2^) is ex- 
plained in greater detail in Appendix |A|. There, we also 
consider the number of qubits, r, that must be used to 
represent the pointer. We show that if the gap is only 
polynomially small in n, it is always sufficient to take 
r = O(logn). However, we argue that generally a single 
qubit will suffice. 



IV. THE GROVER PROBLEM 

The unstructured search problem considered by Grover 
is to find a particular unknown ri-bit string w (the marked 
state, or the winner) using only queries of the form "is 
z the same as u>?" B. In other words, one is trying to 
minimize a function 



h w (z) 



. z = w 

1 , Z ^ W 



(30) 



Since there are 2™ possible values for w, the best pos- 
sible classical algorithm uses 0(2") queries. However, 
Grover's algorithm requires only 0(2™/ 2 ) queries, pro- 
viding a (provably optimal |T^| ) quadratic speedup. In 
Grover's algorithm, the winner is specified by an oracle 
U w with 



U w \z) = (-l) h ^\z) 



(31) 



This oracle is treated as a black box that one can use 
during the computation. One call to this black box is 
considered to be a single query of the oracle. 

In addition to Grover's original algorithm, quadratic 
speedup can also be achieved in a time-independent 
Hamiltonian formulation pQ ] or by adiabatic quantum 
computation (2l], [22]| . In cither of these formulations, the 
winner is specified by an "oracle Hamiltonian" 



H w = l- \w)(w\ 



(32) 



whose ground state is \w) and that treats all orthogo- 
nal states (the non-winners) equivalently. One is pro- 
vided with a black box that implements H Wl where w is 
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FIG. 2: Oracles for the Grover problem, (a) Top: Grover's 
original oracle, (b) Center: An oracle that performs evolution 
according to H w . The double line indicates a classical control 
parameter, the time for which the Hamiltonian is applied, (c) 
Bottom: An oracle that allows one to measure H m . 



unknown, and is asked to find w. Instead of counting 
queries, the efficiency of the algorithm is quantified in 
terms of the total time for which one applies the oracle 
Hamiltonian. 

Here, we show that if we are given a slightly differ- 
ent black box, we can achieve quadratic speedup using 
the measurement algorithm. We let the problem Hamil- 
tonian be Hp — H w and we consider a one-parameter 
family of Hamiltonians H(s) given by @ for some Hb- 
Because we would like to measure this Hamiltonian, it 
is not sufficient to be given a black box that allows one 
to evolve the system according to H w . Instead, we will 
use a black box that evolves the system and a pointer 
according to H w ® p, where p is the momentum of the 
pointer. This oracle is compared to the previous two in 
Fig. H By repeatedly alternating between applying this 
black box and evolving according to Hb (8 p, each for 
small time, we can produce an overall evolution accord- 
ing to the Hamiltonian [sHb + (1 — s)Hp] ®p, and thus 
measure the operator H(s) for any s. 

Now consider the beginning Hamiltonian 



H, 



E 



i - 



(33) 



(i) 

where <7£ is the Pauli x operator acting on the jth 
qubit. This beginning Hamiltonian is a sum of lo- 
cal terms, and has the easy-to-prepare ground state 
l-E'o(O)) = 2~"/ 2 \z), the uniform superposition of all 
possible bit strings in the computational basis. If we con- 
sider the interpolation (||), then one can show || that the 
minimum gap occurs at 



s* = 1 - - + 0{n- 2 ) , 
n 

where the gap takes the value 

ff(s*) = 2 1 -"/ 2 [l + 0(n- 1 )] 



(34) 



(35) 



Naively applying ( p7| ) gives a running time T = 0(2 3n / 2 ), 
which is even worse than the classical algorithm. 

However, since we know the value of s* independent 
of w, we can improve on this approach by making fewer 
measurements. We observe that in the limit of large n, 
the ground state of H(s) is close to the ground state 
\E (0)) of Hb for s < s* and is close to the ground 
state \E (l)) = \w) of Hp for s > s* , switching rapidly 
from one state to the other in the vicinity of s = s*. 
In Appendix |b[ we show that up to terms of order 1/n, 
the ground state and the first excited state \ip~) of 
H(s*) are the equal superpositions 

|^ ± )~-^(|£b(0))±|^(l))) (36) 

of the initial and final ground states (which are nearly 
orthogonal for large n). If we prepare the system in 
the state \Eq(0)) and make a perfect measurement of 
H(s*) followed by a perfect measurement of H(l), we 
find the result w with probability i . The same effect can 
be achieved with an imperfect measurement, even if the 
pointer consists of just a single qubit. First consider the 
measurement of H(s*) in the state |£o(0)). After the 
system and pointer have interacted for a time t accord- 
ing to dio| ) with r = 1, the reduced density matrix of the 
system in the {|^ + ), |VO} basis is approximately 



M''W* cos[ 5 (s*)i/4] 



Jg(s")t/i cos[g(s*)t/4] 

1 



1 

2 V e ~ 

(37) 

If we then measure H(l) (i.e., measure in the compu- 
tational basis), the probability of finding w is approxi- 
mately 



.sin 2 [ 3 ( S *)*/4]. 



(38) 



To get an appreciable probability of finding w, we choose 

t = e(2*/ 2 ). 

This approach is similar to the way one can 
achieve quadratic speedup with the adiabatic algorithm. 
Schrodinger time evolution governed by (0) does not yield 
quadratic speedup. However, because s* is independent 
of w, we can change the Hamiltonian quickly when the 
gap is big and more slowly when the gap is small. Since 
the gap is only of size ~ 2~™/ 2 for a region of width 
~ 2~ n / 2 , the total oracle time with this modified sched- 
ule need only be 0{2 n / 2 ). This has been demonstrated 
explicitly by solving for the optimal schedule using a dif- 
ferent beginning Hamiltonian H' B that is not a sum of 
local terms |2l|, |22| , but it also holds using the beginning 
Hamiltonian (|33|j7 

Note that measuring H(s*) is not the only way to solve 
the Grover problem by measurement. More generally, 
we can start in some u>-independent state, measure the 
operator 



H — H,„ 



K 



(39) 
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where K is also independent of w, and then measure in 
the computational basis. For example, suppose we choose 



K=l-\i>){i>\ 



(40) 



where \ip) is a w-independent state with the property 
~ 2~"/ 2 for all w. (If we are only interested 
in the time for which we use the black box shown in 
Fig. ^](c), i.e., if we are only interested in the oracle query 
complexity, then we need not restrict K to be a sum of 
local terms.) In (40), the coefficient of —I in front of 
\ijj) has been fine-tuned so that + \w) is the ground 
state of H (choosing the phase of \w) so that (w\ ip) is real 
and positive) . If the initial state has a large overlap with 
\ip), then the measurement procedure solves the Grover 
problem. However, the excited state — \w) is also 
an eigenstate of H, with an energy higher by of order 
2-n/2 Thus the time to perform the measurement must 
be f2(2™/ 2 ). 

The measurement procedures described above saturate 
the well-known lower bound on the time required to solve 
the Grover problem. Using an oracle like the one shown 
in Fig. |^(a), Bennett et al. showed that the Grover prob- 
lem cannot be solved on a quantum computer using fewer 
than of order 2™/ 2 oracle queries fLSfl . By a straightfor- 
ward modification of their argument, an equivalent result 
applies using the oracle shown in Fig. ||(c). Thus every 
possible H as in (|3^) that can be measured to find w 
must have a gap between the energies of the relevant 
eigenstates of order 2~™/ 2 or smaller. 



V. DISCUSSION 

We have described a way to solve combinatorial search 
problems on a quantum computer using only a sequence 
of measurements to keep the computer near the ground 
state of a smoothly varying Hamiltonian. The basic prin- 
ciple of this algorithm is similar to quantum computation 
by adiabatic evolution, and the running times of the two 
methods are closely related. Because of this close con- 
nection, many results on adiabatic quantum computa- 
tion can be directly imported to the measurement algo- 
rithm — for example, its similarities and differences with 
classical simulated annealing [p3| . We have also shown 
that the measurement algorithm can achieve quadratic 
speedup for the Grover problem using knowledge of the 
place where the gap is smallest, as in adiabatic quantum 
computation. 

One of the advantages of adiabatic quantum compu- 
tation is its inherent robustness against error p4| . In 
adiabatic computation, the particular path from Hb to 
Hp is unimportant as long as the initial and final Hamil- 
tonians are correct, the path is smoothly varying, and 
the minimum gap along the path is not too small. Ex- 
actly the same considerations apply to the measurement 
algorithm. However, the adiabatic algorithm also en- 
joys robustness against thermal transitions out of the 
ground state: if the temperature of the environment is 



much smaller than the gap, then such transitions are 
suppressed. The measurement algorithm might not pos- 
sess this kind of robustness, since the Hamiltonian of the 
quantum computer during the measurement procedure is 
not simply H(s). 

Although it does not provide a computational advan- 
tage over quantum computation by adiabatic evolution, 
the measurement algorithm is an alternative way to solve 
general combinatorial search problems on a quantum 
computer. The algorithm can be simply understood in 
terms of measurements of a set of operators, without ref- 
erence to unitary time evolution. Nevertheless, we have 
seen that to understand the running time of the algo- 
rithm, it is important to understand the dynamical pro- 
cess by which these measurements are realized. 
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APPENDIX A: THE MEASUREMENT PROCESS 



In Section [II, we discussed the running time of the 
measurement algorithm by examining the measurement 
process. In this Appendix, we present the analysis in 
greater detail. First, we derive the bound on the running 
time by demonstrating (|2|) and (|6|) . We show rigorously 
that these bounds are sufficient as long as the gap is 
only polynomially small and the number of qubits used 
to represent the pointer is r = 0(\og n). Finally, we argue 
that r — I qubit should be sufficient in general. 

Our goal is to find a bound on the final success prob- 
ability of the measurement algorithm. We consider the 
effect of the measurements on the reduced density matrix 
of the system, which can be written as the block matrix 



v X 



(Al) 



where fi = poo, v a = p a0 for a ^ 0, and Xab = Pab 
for a,b ^ 0. Since tr p — 1, /i = 1 — tr \- For ease of 
notation, we suppress j, the index of the iteration, except 
where necessary. The unitary transformation ( |l2| ) may 
also be written as a block matrix. Define e = TS/g. Using 
perturbation theory and the unitarity constraint, we can 
write 



U = 



VU + 0(e 3 ) 
V + 0(e 2 ) 



(A2) 
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where |it| 2 > 1 - e 2 + 0(e 3 ), ||u;|| 2 < e 2 + 0(e 3 ), and V 
is a unitary matrix. We let ||-|| denote the l 2 vector or 
matrix norm as appropriate. Furthermore, let 



1 fct 
k J 



(A3) 



From ([1 
written 



the effect of a single measurement may be 



ft = {UpU r )on. 



(A4) 



where o denotes the element-wise (Hadamard) product. 
If we assume ||^|| = 0(e), we find 



\u\ 2 n - w ] Vv - v ] V^w + 0(e 3 ) 



(A5) 



v = [Vu + ftw - VxV^w + 0(e 2 )] o k . (A6) 

Now we use induction to show that our assumption al- 
ways remains vali d. Initially, = 0. Using the triangle 
inequality in (A6), we find 



IKH < [\W\\+e + 0(e 2 )]k. 



where 



k = max 



ki 3) 



(A7) 



(A8) 



So long as k < 1, we can sum a geometric series, extend- 
ing the limits to go from to oo, to find 



< 



1 



T +0(e 2 ) (A9) 
k 



for all s and for all a / 0. The left hand bound with 
a = 1 gives t > ir/g, which is @). Requiring the right 
hand bound to hold for the largest energy difference gives 
the additional condition 2 r > (E^-i — Eo)/g. Since 
we only consider Hamiltonians H(s) that are sums of 
poly(n) terms of bounded size, the largest possible energy 
difference must be bounded by a polynomial in n. If we 
further suppose that g is only polynomially small, this 
condition is satisfied by taking 



r = O(logra) , 



(A13) 



as claimed at the end of Section |n[ Thus we see that the 
storage requirements for the pointer are rather modest. 

However, the pointer need not comprise even this many 
qubits. Since the goal of the measurement algorithm is 
to keep the system close to its ground state, it would be 
surprising if the energies of highly excited states were rel- 
evant. Suppose we take r = 1; then |k(o;)| 2 = cos 2 (x/2). 
As before, ( p6| ) suffices to make |«oi| 2 sufficiently small. 
However, we must also consider terms involving |Ko a | 2 
for a > 1. The algorithm will fail if the term fiw o k 



in (A6) accumulates to be O(l) over M iterations. This 
will only happen if, for 0(M) iterations, most of ||w|| 
comes from components w a with (E a — Eo)t close to an 
integer multiple of 2ir. In such a special case, changing 
t will avoid the problem. An alternative strategy would 
be to choose t from a random distribution independently 
at each iteration. 



APPENDIX B: EIGENSTATES IN THE GROVER 
PROBLEM 



for all j. In other words, \\i/\\ — 0(e) so long as k is 
bounded below 1 by a constant. 

Finally, we put a bound on the final success probabil- 
ity fS M \ Using the Cauchy-Schwartz inequality in (A5) 
gives 



//>(l-e 2 )/i 



2c- 



0(e 3 



(A10) 



Iterating this bound M times with the initial condition 
fjt(°> = 1, we find 



^ > 1 



M 5 2 



0(Me 3 ) . (All) 



If k is bounded below 1 by a constant (independent of 
n) , we find the condition (|2^) as claimed in Section H. 

The requirement on gives the bound ( p6| ) on the 
measurement time r, and also gives a condition on the 
number of pointer qubits r. To see this, we must inves- 
tigate properties of the function |k(x)| 2 defined in jl5| ) 
and shown in Fig. [l]. It is straightforward to show that 
\k{x)\ 2 < 1/2 for tt/2 < x < n(2 r - 1/2). Thus, if we 
want k to be bounded below 1 by a constant, we require 



tt/2 < [E a (s) - E Q {s)]t/2 < ir{2 r - 1/2) 



(A12) 



Here, we show that the ground state of H(s*) for the 
Grover problem is close to (|3^). Our analysis follows 
Section 4.2 of @. 

Since the Grover problem is invariant under the choice 
of w, we consider the case w = without loss of gener- 
ality. In this case, the problem can be analyzed in terms 
of the total spin operators 



1 



(J) 



(Bl) 



where a = x,y,z and as! is the Pauli a operator acting 
on the jth qubit. The Hamiltonian commutes with S 2 = 



S 2 + S 2 + S 2 , and the initial state has S 2 



(f + l),so 



n ( n 
2 V 2 



we can restrict our attention to the (n + l)-dimcnsional 
subspace of states with this value of S 2 . In this subspace, 
the eigenstates of the total spin operators satisfy 



5 |m a = m) = m\m a = m) 



(B2) 



for 



•m = -f,-f + l,^ 
spin operators and eigenstates, the Hamiltonian is 



^ . Written in terms of the total 



H(s) 



(1- 

+»( 



_ St 

2 



1 - 



n 



(B3) 
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The initial and final ground states are given by |-Eq(0)) 



and \Eo(l)) 



), respectively. 



Projecting the equation H(s)\tj;) 
eigenbasis of S x , we find 



E\tp) onto the 



n 
2 



rp r 



l-sr-A 



(B4) 



where we have defined A = (E — s)/(l — s) and P r = 
2~ n [Z\ . Now focus on the ground state 1-0+) and the first 
excited state \ipJ) of H(s*). By equation (4.39) of ||, 
these states have A± = =F§ 2"™/ 2 (l + 0[l/n)). Putting 
r = in (B4) and taking s = s* from (J3J), we find 



For r ^ 0, we have 

x = | - r ip ; 



V±)(l + 0(l/n)). (B5) 



n v^- 
x(l + 0(l/n)) 



0=1 



Requiring that \ip±) be normalized, we find 



E 



m T = r 

2 



n 

m. = — 

z 2 



-0± 

x(l + 0(l/n)) 



n 

m z = 77 

2 



(2 + 0(l/»)), 



(B7) 



(B8) 
(B9) 



which implies |(m* = f |^±)| 2 = ^ + 0(l/n). From ( |B5|) , 
we also have \{m x — flV'i)! 2 = i +0(l/n). Thus we 
find 



I0d 



in. , 



± 



(BIO) 



(B6) up to terms of order 1/n, which is (|3q). 
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